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ABSTRACT 

We present two-dimensional, integral field spectroscopy covering the rest-frame wavelengths 
of strong optical emission lines in nine sub-mm-luminous galaxies (SMGs) at 2.0<z<2.7. 
The GEMINI-NIFS and VLT-SINFONI imaging spectroscopy allows the mapping of the gas 
morphologies and dynamics within the sources, and we measure an average Ha velocity dis- 
persion of (a) = 220±80kms~ 1 and an average half light radius of (ri/2} = 3.7±0.8kpc. 
The dynamical measure, (V b s /2(T) = 0.9 ±0.1 for the SMGs, is higher than in more quies- 
cent star-forming galaxies at the same redshift, highlighting a difference in the dynamics of the 
two populations. The far-infrared SFRs of the SMGs, measured using Herschel- SPIRlflfar- 
infrared photometry, are on average 370±90M Q yr~ 1 which is ^2 times higher than the ex- 
tinction corrected SFRs of the more quiescent star-forming galaxies. Six of the SMGs in our 
sample show strong evidence for kinematically distinct multiple components with average ve- 
locity offsets of 200±100kms~ 1 and average projected spatial offsets of 8±2kpc, which we 
attribute to systems in the early stages of major mergers. Indeed all SMGs are classified as 
mergers from a kinemetry analysis of the velocity and dispersion field asymmetry. We bring 
together our sample with the seven other SMGs with IFU observations to describe the ionized 
gas morphologies and kinematics in a sample of 16 SMGs. By comparing the velocity and spa- 
tial offsets of the SMG Ha components with sub-halo offsets in the Millennium simulation 
database we infer an average halo mass for SMGs in the range of 13 < log(M[h~ 1 M Q ]) < 14. 
Finally we explore the relationship between the velocity dispersion and star formation inten- 
sity within the SMGs, finding the gas motions are consistent with the Kennicutt-Schmidt law 
and a range of extinction corrections, although might also be driven by the tidal torques from 
merging or even the star formation itself. 



1 INTRODUCTION 



Studies of the Ha nebular emission line have revolutionized our 
understanding of moderate luminosity star forming galaxies at 
z~2. The most extensive study to date has been through the 
Spectroscopic Imaging survey in the Near-infrared with SIN- 
FONI (SINS) of 62 star-form ing galaxies at 1.3 < z < 2.6 
dForste r-Schreib er et al.l 12009b . a subset of known Ha emitters 
from the UV-con tinuum selected galaxies dSteidel et al.l |2004 
lDaddietal.ll2004h . In the SINS galaxies, the Ha velocity fields 
are mapped to trace the dynamics of the ionized gas within 
the galaxies, with a striking result that the star-forming galaxies 



* Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important partic- 
ipation from NASA. 



can be split into three equally-populated groups; turbulent disks 
which are rotation dominated, compact sources which are disper- 
sion dominated, and merging systems. These galaxies populate 
the 'main sequence' of star-form ing galaxies at these redshifts (eg 
iDaddi et al.|[2007l : lNoeske et alj|2007n . The SINS survey, however, 
only has 6% of galaxies with inferred star-formation rates (SFRs) 
greater than 5OOM0yr _1 , which is the domain of galaxies se- 
lected through bright far-infrared emission redshifted into the sub- 
millimetre wavelengths, or sub-millimetre galaxies (SMGs). Fur- 
ther, the poorly constrained dust extinction corrections in the SINS 
galaxies implies that this 'high-SFR' tail may actually have sub- 
stantially lower SFRs. There is thus an obvious need to explore the 
nebular emission line dynamics of the SMG population to compare 
to the 'main sequence' of star-forming galaxies. 

Since their discovery, SMGs have been the focus of intense 
multi- wavelength study. With 15 years since the first sub-mm sur- 
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veys (eg ISmail et al.|[l997l ; Hjughes et alj|l998h our understanding 
of this population has grown considerably, although still limited 
by only ha ndfuls of objects studied in detail at hig h spatial res- 
olutio n (eg ISwinbank etall 12004 lEngel et alj l201fA ICarilli et al.l 
l20ld) . These dus ty galaxies have a r e dshift distribution w hich 
peaks at z ~2.5 dChapman et ai]|2005l ; IWardlow et all 1201 lh . At 
these redshifts, the observed sub-mm flux densities imply far- 
infrared luminosities >10 12 L0 (Ultra-Luminous Infra-Red Galax- 
ies - ULIRGs). With space densities at z ~2 a factor ~1000x 
higher than today dChapman et al.l 120031) . these luminous sys- 
tems contribute ~30% of the star-formation activity at z ~2-3, 
compared to <1% locally. It is noteworthy that some confirmed 
ULIRGs at z ~2 have hotter SED shapes resulting in 850 pirn 
fluxes lying below the sensitivities of current surveys, and the high- 
redshift ULIRG population co uld be up to twice as large as im- 
plied by 850 nm surveys alone dChapman et al.l2004bl ; ICasev et al.l 
l2009l ; lMagdis et al]|2010h . These galaxies which have been missed 
by 85 fxm surveys, have been detected by BLAST jCasev et al.l 
l201lh and Herschel dMagnelli et afll20ld : IChapman et alj|2010h , 
and confirmed as 'SMGs' with 250-500 /im flux densities still 
implying ULIRG-like luminosities. The intense star-formation 
rates (SFRs ^lOOOMn yr -1 ) and large molecular gas reservoirs 
(~10 n Mir dGreve et alj2005t Bothwell et al. in prep.) inferred for 
the SMGs means that they could build a stellar mass of 10 11 M in 
just 100 Myr. The space densities and inferred duty cycles of SMGs 
suggest they r epresent the formation phase of today's massive ellip- 
tical g alaxies dLillv et al.ll999l ; lGenzel et all2003HSwinbank etall 
12008) . 

SMGs are clearly an influential population of early galaxy 
evolution, but attempts to model their basic properties have proved 
problematic. The most successful semi-analytic models have had 
to radically alter their prescriptions for starbursts during merg- 
ers in order to account for SMGs whilst r emaining self-consi sten t 
with other constraints dBaugh et al.ll2005h . lNaravanan et"al1 d2009l) 
looked at isolated conditions in hydrodynamical simulations that 
could lead to very high SFRs, and proposed a me rger-driven for- 
mation scenario for SMGs. Subsequent analysis bv lHavward et al.l 
(2011b) showed that these merging galaxies may represent the 
bright end of the SMG number counts, but found that the fainter tail 
of SMGs (Sgso ~3 mjy) are t ypical ly made up of isolated galaxies 
and galaxy pairs. iDave et al.l d2010l) identified SMGs in their cos- 
mological simulations as typically very massive galaxies which are 
being harassed by small fragments, offering an alternative possi- 
bility to the major-merger picture for SMGs. Observations of indi- 
vidual SMGs at high-resolution are therefore required to guide the 
models and constrain the driver for the high SFRs. 

High-resolution radio and CO observations have shown that 
the star-formation and gas can be extended on scales up to ~ 10 kpc 
in SMGs, and multiple components are common within these 
extended distributions dCh apman et al.l l2004d: iBiggs & IvisorJ 
20081: iTacconi et al.ll2008tlCasev et alJl2009l ; lBothwell et alJl2o7q 
Ivison et al.l 1201 lh . These sizes are much larger than the typi- 



cal ~lkpc extent of local ULIRGs (eg ISakamoto et alj 120081) 
whose intense SFRs are of t en triggered by mer gers and interactions 



vnose intense i^Ks are otten triggered by mergers and interactions 
Sanders & Mirabel| [l996). Eng el et al.l feoiOh and lBothwell et alj 



2010) concluded that most SMGs are formed in major mergers 



many exhibiting distinct kinematic components and disturbed gas 
morphologies (indicative of pre-coalesced to late stage mergers) 
from their high resolution interferometric CO line studies . 

Studying SMGs in Ha complements the CO studies, as lower 
levels of star formation are detectable through Ha emission than 
through the existing CO or radio measurements (from the IRAM- 



Plateau de Bure Interferometer and the Very Large Array respec- 
tively). Further, merging galaxies may not induce similarly high 
levels of star formation in all components of the merger, and Ha 
may elucidate a more exten sive merging system than apparent in 
CO. ISwinbank et alj d2006l) find multiple components in four of 
the six SMGs they studied using integral field unit (IFU) spec- 
troscopy to map Ha, with average velocity offsets of ~180kms~ 
and spatial offsets up to ~8kpc. They concluded that SMGs are 
similar to the local, more compact ULIRGs in that they exhibit 
multiple components. However, they were not able to resolve kine- 
matic information on scales <5kpc (i.e., within each component) 
and thus the gas dynamics of the multi-component systems could 
not b e established. Using higher spati al resolution IFU observa- 
tions, |MenendezjDelmestte|et^] |20l2) measure the gas kinemat- 
ics within the multiple components of three SMGs. However the 
majority of the Ha fields of these three sources are dominated by 
broad-line AGN, and only small sections of the narrow line regions 
have [NII]/Ha ratios which are consistent with photo-ionization 
from star formation. There has thus remained a need to further ex- 
plore a larger sample of SMGs with high resolution, sensitive IFU 
spectroscopy, including the large fraction of SMGs (>50%) which 
do not show such obvious indications of strong AGN in their opti- 
cal spectra. 

In this paper we present the gas dynamics of nine SMGs ob- 
served with the Gemini-NIFS and VLT-SINFONI. In Section[2]we 
present the sample and the observations. In Section [3] we analyse 
the spatially integrated and resolved observations. The dynamical 
properties of the SMGs are presented in Section |4] along with the 
comparison to other galaxy surveys and the kinemetry analysis. In 
Section[5]we discuss the merger properties and use the component 
properties to infer the average halo mass. Finally, in Section [6] we 
give our main conclusions. All calculations assume a flat, ACDM 
cosmology with = 0.7 and Ho = 71kms~ Mpc -1 in which 
0.5" corresponds to a physical scale of ~4 kpc at z — 2 (the median 
seeing and redshift of our observations and targets respectively). 



2 SAMPLE SELECTION AND OBSERVATIONS 
2.1 Sample 

Our IFU survey comprises nine radio-identified ULIRGs (listed 
in Table [T}, targeted from various surveys, preferentially chosen 
from our equatorial and southern fields to enable further study 
with ALMA. //erac/ie/-SPIRE fluxes are available for the major- 
ity of the sources (Table Q3 enabling the far-infrared luminosities 
of the samples to be derived directly. We find that all sources 
have Lfir >1O 12 L . The targets designated 'SMM' are SMGs 
with Sgso/jm > 4mJy, while those listed as 'RG' are unde- 
tected at 850 /im, with typical 3a rms of 4.5 mjy. These sources 
are, however, still found to be bright at 250, 350, & 500 /im 
and therefore represent the ULIRGs with hotter dust tempera- 
tures causing the 850 /im fluxes to fall below the detection limits. 
Eight of the sources have spectroscopic redshifts derived from our 
longslit detections of the Ha emission line. For the ninth source, 
SMMJ2217+0017, we obtained the redshift directly with our IFU 
observations. SMMJ2217+0017 was identifie d as a bright SMG 
(Snoi^m = 4.9 ± 0.7 mj y) in the SA22 field dTamura et al.l2009l) 
with a VLA counterpart dChapman et alj l2004ah . and an optical 
through mid-IR (Spitzer-IKAC) photometric redshift constrained to 
2 ~2-2.5. 

The SMGs span a redshift range of 2, 0< z <2.7, which cov - 
ers the median of the parent sample from IChapman et ail d2005l) . 
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Figure 1. The H« luminosities from our IFU observations (not cor- 
rected for extinction) versus radio luminosities of the sour ces studied in 
this w ork, compared to the sample of SMGs studied in Swinba nk et alj 
l2004)using long-slit spectroscopy. The SFRi ,4GHz ax is represents the 
SFRs derived from the rad io luminosities, using the far-infrared-r adio cor- 
relation Jlvis on et al. 2010) and the conversion to SFR from Kennicutt 



1998a). The SFRh q axis assumes the conversion to SFR from 



Kennicutt 



1998a) and both SFRi,4ghz and SFRhq conversions include a factor 



of 1.7 to convert from the Salpeter to Chabrie r IMF. The sample of SMGs 
with p reviously published IFU o bservations fewi nbank et al. 2005), l2006r . 
Menendez-Delmestre et al. 2012) are highlighted. We also plot the Hot lu- 
minosity for RGJ0332-2758 ICasev etai]|201lh which has IFU observa- 
tions presented in Casey et al. (in prep.). Our targ et sample has a simila r 
range of radio and Ha luminosities to the SMGs in Swinban k et alj J2004I) . 
The dotted line denotes equality in the two SFR indicators whereas the 
dashed line correspond to the relationship between the two SFR indica- 
tors for an extinction o f A v = 2.9 ± 0.5 (the average value for SMGs 
from iTakata et al.l l 2006 ) which decreases the observed SFR^f^ relative to 
the radio. We also mark the target sources which potentially host AGN, as 
discussed in Section 3.3 and the previously published sources which are 
known to host AGN. 



and places Ha in the if-band atmospheric window. Fig. Q] shows 
the range of 1.4 GHz radio (referenced in Table [T} and Ha lu- 
minosities in the targ et sample compared to th e larger sample 
of SM Gs studied in ISwinbank et alj d2004h and IChapman et alj 
J2005h . demonstrating our IFU sample is representative. We use 
the radio luminosities and radio-derived SFRs in Fig.Q] rather than 
the directly measured far-infrared luminosities and SFRpir val- 
ues, to enable a comparison to the complete parent sample. Fig. Q] 
also shows that our sample is consistent with having extinctions of 
A„ ~2.9, corresponding to the a verage A„ measure d for a sample 
of SMGs from the Ha/H/3 ratio (Takata et al. 2006), which brings 
the Ha SFRs in line with those inferred from the far-infrared. 



it into 29 slices of width 0.103". The dispersed spectra from the 
slices are reformatted on the detector to provide 2-dimensional 
spectro-imaging, in our case using the 7f-band grism covering 
a wavelength range of 2.00-2.43 /im. The observations were 
performed using an ABBA configuration in which we chopped by 
10" to blank sky to achieve sky subtraction. Individual exposures 
were 600 s and each observing block was 2.4 ks, which was 
repeated between three and four times resulting in the on source 
integration times (which do not include the sky observations) given 
in TableCD 

SMMJ2217+0017, SMMJ0217-0503, RGJ0332-2732 and 
SMMJ0217— 0459 were observed with the Spectrograph for IN- 
tegral Field Observations in the Near Infrared (SINFONI) on the 
VLT between 2011 July 13th and 2011 September 23rd as part of 
programme 087.A-0660(A). All of our observations targeted the 
redshifted Ha and [Nll]AA6548. 1,6583.0 emission lines in the K- 
band. We use the H+K grating on SINFONI which gives a spectral 
resolution R=1500 covering a range of 1.45—2.45 fim. The field of 
view is sliced into 32 slices each of which is imaged onto 64 pix- 
els. We chose the width of the slices to be 0.25" to provide a 8" x 8" 
field (obtaining 32 by 64 spectra of the image with each pixel cov- 
ering 0.25" by 0.125"). We move the target around 4 quadrants of 
the IFU for sky subtraction purposes, however retaining the target 
in the field of view, therefore the effective field of view is 4"x4" 
forSMMJ2217+0017, RGJ0332-2732 and SMMJ0217-0459 and 
4" x 8" for SMMJ0217-0503 where we only use 2 quadrants in or- 
der to observe the entire extended system. The individual exposure 
times are 600s with resulting on source exposure times (which do 
not include the sky observations) given in Table [T] 

We reduced the NIFS data with the standard Gemini IRAF 
NIFS pipeline and the SINFONI data with the SINFONI ESOREX 
pipelines both of which include extraction, sky subtraction, wave- 
length calibration, and flat-fielding. To accurately align and mo- 
saic the individual datacubes we created white light (wavelength 
collapsed) images around the redshifted Ha emission line from 
each observing block and centroid the galaxy within the data cube. 
These were then spatially aligned and co-added using an aver- 
age with a 3ct clipping threshold to remove remaining cosmetic 
defects and cosmic rays. Flux calibration was carried out by ob- 
serving bright A0V standard stars at similar airmass to the target 
galaxy immediately after each observing block for both the NIFS 
and SINFONI observations, and we then also compare the inte- 
grated fluxes to those measured from longslit observations of the 
targets. From the reduced standard star cubes, we measure a see- 
ing of FWHM~0.5" for both the NIFS and SINFONI observations 
which corresponds to ~4kpc at z ~ 2. The spectral resolution of 
the data (measured from the sky-lines at ~2.2/^m) is cr=1.8 A or 
25 kms -1 for the NIFS observations and <j=6.3 A or 88 kms -1 for 
the SINFONI observations. In all following sections, quoted line 
widths are corrected to take account of the instrumental resolution, 

Ocorr = yip^bs ~ a sky)- 



2.2 Observations and Reduction 

Observations of RGJ0331-2739, SMMJ0217-0505, 
SMMJ0333-2745, SMMJ2218+0021 and SMMJ0332-2755 
were made with the Gemini-North/Near-Infrared Integral Field 
Spectrograph (NIFS) between 2010 September 24th and 2011 
February 28th as part of program GN/2010B/42. The Gemini-NIFS 
IFU uses an image sheer to take a 3.0" x 3.0" field and divides 



3 ANALYSIS 

3.1 Spatially integrated spectra 

In Fig. [2] we show the integrated one-dimensional spectra of the tar- 
gets with strong Ha detections from our IFU observations. We do 
not include the weak Ha detection in SMMJ0217— 0459 since the 
IFU observations are low S/N and do not provide further informa- 
tion than the previous longslit observations. We use the integrated 



© 0000 RAS, MNRAS 000. 000-000 



4 S. Alaghband-Zadeh et al. 



IS ourcc 


R A 

tSJ\ 


UCL 


ZHa 


il.4GHz 


•J250nm 


"J350fim 


i>500fim 




tint 


IN U Lc-S 




(J2000) 


(J2000) 




Ouiy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(ks) 




RGJ0331-2739 1 


03:31:52.82 


-27:39:31.5 


2.3429[1] 


965 ± 16 


36.0 ±3.6 


47.0 ±4.2 


42.0 ±4.8 


1.9 ± 1.5 


4.8 


1 


SMMJ0217-0505 


02:17:38.92 


-05:05:23.7 


2.5305[2] 


41 ± 11 


24.8 ±3.6 


26.2 ±4.2 


26.0 ±4.8 


7.1 ± 1.5 


4.8 


2 


SMMJ0333-2745 


03:33:15.43 


-27:45:24.4 


2.6937[2] 


75.8 ±6.9 


31.2 ± 3.6 


41.0 ±4.2 


35.3 ±4.8 


9.2 ±1.2 


4.8 


3 


SMMJ22 18+0021 


22:18:04.42 


00:21:54.4 


2.5171[4] 


44 ± 10 








9.0 ±2.3 


4.2 


4 


SMMJ0332-2755 


03:32:43.20 


-27:55:14.5 


2.1236[3] 


106 ±9 


32.8 ±3.6 


31.6 ±4.2 


21.2 ±4.8 


5.2 ± 1.4 


3.6 
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SMMJ2217+0017 1 


22:17:42.25 


00:17:02.0 


2.2777[2] 


58 ±11 










9.6 
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RGJ0332-2732 1 


03:32:56.75 


-27:32:06.3 


2.3152[3] 


160 ± 20 


27.5 ± 3.6 


27.5 ±4.2 


20.3 ±4.8 


2.4 ± 1.4 


2.4 


7 


SMMJ0217-0503 2 


02:17:38.62 


-05:03:39.5 


2.0280[3] 


185 ± 12 


36.3 ±3.6 


33.9 ±4.2 


21.7 ±4.8 


4.4 ± 1.7 


4.8 


2 


SMMJ0217-0459 


02:17:25.12 


-04:59:37.4 


2.3090[5] 


57 ± 10 


21.7 ±3.6 


22.2 ±4.2 


17.0 ±4.8 


4.5 ± 1.9 


7.2 
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Table 1. Table of the SMG fluxes, positions, redshifts and on source exposure times. The redshifts quoted are derived from the IFU spectra measurements of 
the Ha emission line. The values in [] in the redshift column re pre sent the error on th e la st decimal place . The r adi o flux and position r efer ences are given 
in the fi nal column. 1 : llyison et al J Boidb , 2: llvison et all l2007t) , 3: lBiggs et aflfeOllh . 4: IChapman et alj 120051) . 5: lDunlop et a?] J2010h . 6: IChapman et alj 
l2004j) . 7: lMiller et alj i2008l) . 1 A GN candidate, and 2 broad line AGN, a s discussed in Section 3.3. The S250nm, S3S0iJ.m and Ssoo^m fluxes are measured 
from Herschel-SPWE observations I Griffi n et alj2010l:IOliver et al.l2012l) . The errors quoted are the source extraction errors but we n ote these are dominate d 
by the confusion no i se (5.8, 6.3 and 6.8 mjy/beam for S250;zm, S35U(zm and Ssoo/jm respectively). The Ssso^m fluxes are taken from Ch apman et afl i2005): 
ICoppin et al] l2006h : IWeiB et alj J2009I) . 



spectra to measure basic properties of the sample (z, Sua and <th q ) 
detailed in Table [2] We fit Gaussian profiles to the Ha and [Nil] 
emission lines simultaneously in the integrated spectra. We also at- 
tempt to fit a broad Ha component however we find in all cases that 
the there is no improvement to the \ 2 °f this fit compared to the fit 
without the broad component so we do not use the broad compo- 
nent fit. Errors on these values were calculated by altering the fit, 
changing one parameter at a time and allowing the other parameters 
to minimise, until the \ 2 changed by 1 . We only use the flux from 
the pixels in which Ha is detected at >3cr to create the integrated 
spectra. For the faint SMMJ0217-0459 we quote the Ha redshift 
(which confirms the redshift gained from previous Ha slit observa- 
tions), Ha flux and a but we are unable to measure any resolved 
properties. For SMMJ0217— 0503b we fit a double Gaussian to the 
Ha profile since it shows a double-peaked disk-like profile. Indeed 
the disk-like profile of SMMJ0217— 0503b is confirmed in Section 

EU 

We measure an average Ha flux of < Sh q >= 3 ± 1 x 
10~ 16 ergs -1 cm~ 2 and an average velocity dispersion of (a) = 
220 ±80 kms-Vlnall sources we clearly detect [NII]A6583 emis- 
sion except SMMJ0332— 2755 for which we quote a 3u upper limit 
to the ratio of [NII]/Ha. For our total sample, we measure an av- 
erage [NII]/Ha ratio of 0.52±0.07 in the sources where [Nil] was 
detected. The large spectral coverage of SINFONI allows for mea- 
surements of H/3 in the three SINFONI sources. H/3 however lies 
in a noise dominated part of the spectra such that no useful limit 
can be placed on the ratio of Ha/H/3. We therefore use an extinc- 
ti on value of A v = 2. 9 ± 0.5, taken as the average extinction given 
in lTakata et alj {2006), to correct the Ha luminosities since these 
are dusty systems. The correcte d luminosities we re converted to 
SFRs using the conversion from lKenmcuttl (l998b) which includes 
a con version factor of 1.7 t o adjust from a Salpeter IMF (|Salpeter 
Il955l) to a Chabrier IMF dChabrienl2003l). We app l y the extinc- 
tion correction as detailed in Forster-Schreiber et al. I (l2009l) . using 
the correction at the wavelength of Ha of Ahq = 0.82A V . The 
total correction applied is therefore io°- 33Av . We provide the un- 
corrected Ha fluxes in Table[2]along with the corrected SFRs. 

Since we do not have individual extinction estimates to the 
Ha line, we derive far-infrared luminosities of our SMGs in or- 
der to provide unextincted estimates of the SFR. For seven of our 
SMGs we can use the //er.sc/;e/-SPIRE observations from the DR1 



data release dGriffin et ail 1201 fi lOliver et al.1l2012l) to derive the 
far-infrared luminosities. We extract the SPIRE fluxes at the po- 
sitions of our SMGs (Table [TJ and fit modified black-bodies to 
the photometry at the known redshift. In this analysis, we fix the 
dust emissivity, /3 = 2 (Maanelli e t al.l2012h . but allow the dust tem- 
perature and bolometric luminosity to vary. For SMMJ22 18+0021 
and SMMJ2217+0017, for which do not have 250, 350, & 500 jum 
fluxes, we scale the far-infrare d SEP of SMM J2 135-0 102 (the well 
studie d lensed SMG at z=2.3: ISwinbank et alj|2010bl : llvison et alj 
1201 lh . to the 850 /im flux for SMMJ22 18+0021 and the 1.1mm 
flux for SMMJ2217+0017 (Si.imm = 4.9 ± 0.7 mJy). We com- 
pare the radio luminosities to the far-infrared luminosities to probe 
the far-infrared ra dio correlation. We derive, q using Equation Q] 
dHelouet all 19851) . 

, / Lfir \ /Li.4GHz\ 

q = l0g v 9.8x10-i 5Lq ) l0g { WHFtJ (1) 

We find an median q of 2.0±0.2, for the sources without AGN 
signatures, and a similar q=1.9±0.5 for our entire sample, which 
are consisten t with the q value of a sample of SMGs studied in 
Kovacs et al] 12006J) (q~2.1). We convert the derived far-infrared 
luminosities fro m the SEP fits to SFRs using Equation[2] which is 
the relation from Ken nicutJ Jl998al) adjusted from a Salpeter to a 
Chabrier IMF using a factor of 1.7. 

SFR(M yr" 1 ) = Lpmtergs" 1 ) x 2.6 x 10~ 44 (2) 

The SFRfir values are, on average, approximately in line with 
the extinction corrected SFRn a values indicating that using A v = 
2.9 ± 0.5 brings the two SFR indicators into agreement. 

3.2 Spatially resolved spectra 

We use the spatially resolved Ha emission line maps from our 
NIFS and SINFONI IFU observations to estimate the distribution of 
star-formation and dynamics within the galaxies. After constructing 
the datacubes, we proceed to fit the Ha emission line in each pixel. 
The spectra were averaged over 3x3 spatial pixels (0.15" x 0.15" 
for the NIFS data and 0.375" x 0.375" for SINFONI), increas- 
ing this up to 5x5 pixels (0.25" x 0.25" for the NIFS data and 
0.625" x 0.625" for SINFONI) if the signal was too low to give a 
sufficiently high S/N. In regions where this averaging process still 
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Figure 2. ID integrated spectra of sources showing the location of the Ha line at 6562.8 fim (blue dotted line) and the [Nil] lines at 6583.0,6548.1 (red 
dotted lines). The NIFS spectra have been smoothed by 2.3 pixels to match the resolution of the SINFONI spectra. All spectra are continuum subtracted. The 
dashed lines represent the fit to the Ha and [Nil] lines. The large ratios of [NII]/Ha in SMMJ2217+0017 and RGJ0332-2732 lead us to suggest that these 
sources may host AGN. We mark RGJ033 1—2739 as an AGN candidate' due to its high radio luminosity, however the low value of [NII]/Ho and further 
analysis of the Ha velocity field suggest the Ha properties are not being driven by the AGN activity. SMMJ0217— 0503 is divided into its two clear components 
(as can be seen in the maps of Ha). SMMJ0217— 0503a is most likely an AGN with broad Ha emission blending with the [Nil] lines. SMMJ0217— 0503b has 
a disk-like, double-peaked Ha profile. 



Source 


Shq ""Ha 


SFR Ha 


SFRi. 4G Hz 


SFR FIR 


[NIIJ/Ha 


v obs 


ri/2 


Kasym 




(10- 16 ergs" 1 cm- 2 ) (kms^ 1 ) 


(Mq yr- 1 ) 


(Moyr- 1 ) 




(kmr 1 ) 




(kpc) 





RGJ0331-2739 1 


3.7±0.3 


83±8 


600±200 


8000±100 


280±60 


0.41 ±0.08 


120±20 


2.99±0.04 


0.60±0.05 


SMMJ0217-0505 


2.1±0.2 


130±2O 


500±200 


400±100 


210±50 


0.6±0.1 


400±50 


2.61±0.05 


0.551 ±0.007 


SMMJ0333-2745 


2.1±0.2 


70±10 


500±200 


870±8O 


340±50 


0.3±0.1 


150±30 


2.05±0.06 


5.5±0.3 


SMMJ22 18+0021 


1.8±0.3 


170±30 


400±200 


400±100 


500±100 


0.2±0.1 


380±50 


2.88±0.06 


1.7±0.1 


SMMJ0332-2755 


2.8±0.8 


200±100 


400±200 


610±60 


19O±50 


<0.7 


360±40 


2.48±0.1 


lldbl 


SMMJ2217+0017 1 


0.75±0.06 


150±10 


130±50 


450±90 


900±100 


0.7±0.1 


300±100 


4.8±0.6 


0.61±0.05 


RGJ0332-2732 1 


1.7±0.2 


180±2O 


300±100 


1300±200 


200±40 


0.8±0.1 


200±100 


3.1±0.3 


1.3±0.2 


SMMJ0217-0503 


9.2±0.5 


740±40 


1200±400 


1110±70 


190±50 




800±100 


9.0±0.3 


3.9±0.3 


SMMJ0217-0503a 2 


6.4±0.3 


740±40 


800±300 








300±100 


4.2±0.3 




SMMJ0217-0503b 


1.0±0.1 


90±20 


120±50 






0.4±0.1 


600±100 


4.6±0.3 


0.249±0.005 


SMMJ0217-0459 


1.3±0.4 


100±100 


200±100 


460±80 


150±50 


0.4±0.3 









Table 2. Table of the SMG properties measured and derived from the integrated spectra and maps. The SFRhc values are calculated using Kennicutt 1 1998a) 
and include a conversion factor of 1. 7 to adjust to a C habrier IMF. The SFRhc values are extinction corrected using an average extinction of dusty, star 
forming ULIRGs o f A v = 2.9 ± 0.5 <Takata et al.l2006l) . The SFRi .4ghz values are calculated from the radio fluxes using the far infrared-radio correlation 
llvisonetal]|2010l) . The SFRpir values are calculated from the SED-derived far-infrared luminosities. V b s is the difference between the minimum and 
maximum velocities measured across the Ha velocity field (taken as the 5th and 95th percentiles of the velocity distribution). We quote the properties of the 
whole SMMJ0217— 0503 system as well as the two clear components. The properties of SMMJ0217— 0503b assume a double Gaussian fit to the disk-like Ha 
emission line profile. The broad line fit to Ha of SMMJ0217-0503 and SMMJ0217-0503a (the AGN component) covers Ha and the blended [Nil] doublet, 
therefore it is not possible to measure the [Nil] flux. Caution should be taken with the SFRi 4GHz value for RGJ033 1—2739 since the high radio luminosity 
suggests the radio flux is dominated by the AGN and therefore gives an unphysically high SFR using the far-infrared-radio correlation. The ratios of [NII]/Ha 
we quote include the flux fro m the 6583/xm [Nil ] line only since only this line is included when classing systems as being potentially AGN dominated using 
the comparison to the data in lKewlev et al. 12006). J AGN candidate and 2 AGN as discussed in Section 3.3. 
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failed to give an adequate S/N, no fit was made. Using a contin- 
uum fit, we required a minimum S/N of three to detect the lines, 
and when this criterion is met, we fit the Ha emission line with 
a Gaussian profile allowing the normalization, central wavelength 
and width to vary. 

We simultaneously fit a double Gaussian to the [Nil] lines with 
the fit to Ha. If the improvement in the \ 2 of this fit, compared to 
the single Gaussian fit to Ha, is such that A\ 2 > 9 (corresponding 
to 3a) then the Ha + [Nil] doublet fit is used instead. To calculate 
the error in the line parameters, we alter the fit parameters one at a 
time and allowing the other parameters to minimise until A\ 2 =l, 
corresponding to a formal la error. Figs.|3]and|4]show the results 
from the linefitting; 2-dimensional maps showing the Ha intensity, 
velocity and velocity dispersion fields. The yellow circles represent 
the FWHM of the standard stars showing an estimate of the see- 
ing dominated PSF. This estimate is from shorter exposures (~10s) 
than the science exposures and therefore underestimates a possi- 
ble PSF broadening in the long exposure science observations. We 
use the Ha intensity maps to measure the half-light radii, ri/2, by 
summing the Ha flux within increasing circular apertures, centred 
on the middle pixel of the Ha intensity distribution, until the flux 
enclosed is equal to half the total flux. The ri/2 values are then 
corrected for the PSF. 

Figs. [3] and [4] show that in all sources the Ha distribution is 
extended on scales >5 kpc with an average Ha half-light radius of 
( r i/2) — 3.7 ± 0.8kpc. We com pare our Ha ri/ 2 me asurements 
to molecular gas sizes in SMGs. iTacconi et alj |2006) find a me- 
dian ri/2 of <2±0.8kpcQ in their sample of SMGs observed in 
CO which i s in ag reement with additional CO measurements from 
lEngel et al] d2010l) (median equivalent radius of 2. 2± 1.4 kpc) de- 
spite the fact that these measurements cover a range of high-J CO 
transitions possibly tracing different distributions of the gas. These 
CO sizes are also comparable to the 1 .4 GHz radio siz e s for SMGs 
found by Chapm an et alj d2004ch and|giggs & Iviso (me- 
dian equivalent radius of 1.6±0.4kpc). These compact sizes are 
lower than the average half-light radius we measure in the Ha sam- 
ple (3.7±0.8kpc) suggesting that the Ha observations trace more 
of the gas distribution in the SMGs. However we do not have CO 
line data and Ha observations in the same sample of SMGs there- 
fore comparing the two samples can only be done statistically. 

All sources, except SMMJ2217+0017 and RGJ0332-2732, 
exhibit multiple peaks in the Ha intensity distributions with two 
sources showing clear separations between multiple components 
(SMMJ0332-2755 and SMMJ0217-0503). The velocity and ve- 
locity dispersion fields are not smooth and do not show evidence 
for disk-like profiles which is further addressed in Section |4~2l 

3.3 AGN identification 

Previous IFU studies of SMGs have ma inly targeted galaxies 
show i ng strong AGN in their Ha spec tra dSwinbank et alj 120051 
120061 : iMenende z-Delmestre et alj|2012l) . whereas the majority of 
our sample do not display clear Ha AGN signatures. Indeed, the 
sub-mm and/or far-infrared detections of all the sources in our 
sample provide a baseline of evidence that they are dominated by 
star formation (implied ~1000M Q /yr). In SMGs, AGN are found 
to only contribut e at low levels (<20%) to the bolometric output 
dAlexandej|2004 . However, even if an AGN is not dominating the 

1 All given measurements of FWHM and lineal' diameter are converted to 
equivalent half-light radius by dividing the diameters by 2. 



bolometric output it may still affect the Ha distribution and dynam- 
ics, being energetic e nough to drive ionized gas over scales of a few 
kpc (Nesvadba et al. 2008; Harrison et al. in prep.). 

Two of our sources (RGJ0332-2732 and SMMJ2217+0017) 
have high ratios of [NII]/Ha in their integrated spectra (Fig. 
O, indicating that these sources may host an AGN (a limit of 
[NII]/ Ha >0.7 is typically adopted for AGN dominated systems 
dKewlev et al]|200d) ). Given the high SFRs, it is also possible that 
the high ratios are caused by shock ionization from a galactic wind, 
as oppose d to an AGN ionization fie ld, raising the observed value of 
[NII]/Ha dvan Dokkum et aljfeooll) . The detection of Ha extended 
over a large scales (~ 10 kpc), with no obvious outflow kinematics, 
further suggests these systems may still be star formation domi- 
nated in their Ha kinematics. 

We expect the high [NII]/Ha ratio from AGN excitation would 
typically be concentrated to the inner few kpc. We search the off- 
nuclear regions of RGJ0332-2732 and SMMJ22 17+00 17 for gra- 
dients in the [NII]/Ha excitation, finding that SMMJ22 17+00 17 
shows very concentrated [Nil] ([NII]/Ha=0.7±0.1), with low 
[NII]/Ha ([NH]/Ha=0.3±0.1) in regions around the nucleus, likely 
indicating a compact AGN core surrounded by star formation. 

RGJ0332— 2732, by contrast, shows continued strong 
[NII]/Ha outside the central nucleus and observations of the [OIII], 
H/3 emission indicate photoionization by an AGN also (Harrison et 
al. in prep.). The ratio of [NII]/Ha is higher in the northern region 
of this source ([NII]/Ha=0.9±0. 1 rather than [NII]/Ha=0.7±0. 1 
in the southern region). It is therefore possible that the AGN may 
be ionising the gas to the north only and the remaining extended 
Ha dynamics we observe are unaffected by AGN activity. The 
[Nil] velocity field follows the same galaxy-wide rotation pattern 
as the Ha. However, since there is evidence for the velocity fields 
in RGJ0332-2732 and SMMJ2217+0017 to represent gas which is 
photoionized by an AGN, given the high [NII]/Ha values, we mark 
these sources as an AGN candidates' in the analysis. 

The high radio flux of RGJ033 1-2739 
(Si.4GHz=965±16 /iJy) suggests the presence of an AGN. 
Although it is luminous in the far-infrared as well, the radio 
luminosi ty exceeds that expected from the far-infrared-radio 
relation dlvisonetalfeoiOt) by a factor 29. We do not see evidence 
for any AGN in the Ha, [Nil] fields, and the AGN-radio core could 
be very extincted, buried in the system. It is clear from the velocity 
fields of RGJ033 1-2739 (Fig. [3] and Section O that merger 
kinematics appear to dominate over any possible AGN-driven 
properties, and the [NII]/Ha ratio lies co mfortably withi n that 
expected for a z ~ 2 star forming galaxy dErbetal.1120061) . We 
therefore do not consider this source as an AGN in the remaining 
analysis for purposes of analyzing its nebular line properties. 

SMMJ0217— 0503 is a particularly complex system. It is clear 
in Fig.|4]that SMMJ0217-0503 consists of two spatially separated 
merging galaxies, with integrated spectra for the individual compo- 
nents shown in Fig. [2] The northern source (SMMJ02 17— 0503a) 
reveals a broad Ha line indicating a dominant AGN. The south- 
ern source (SMMJ0217— 0503b) appears to be star formation dom- 
inated (low [NII]/Ha) with a disk-like velocity field and a double 
peaked Ha line profile in the integrated spectrum (Fig|2}. How- 
ever, both components in this~20kpc-separated major merger ap- 
pear to be far-infrared-luminous (e.g., Bothwell et al. in prep.; 
Alaghband-Zadeh et al. in prep). Although we measure the proper- 
ties of these two sources in SMMJ0217— 0503 separately in the fol- 
lowing analyses, we consider the whole system as a major merger. 
We note however that the large velocity offset of the two compo- 
nents (670±70kms~ 1 ) implies that it is possible that this is a flyby 
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Figure 3. Maps of the properties of the Ha emission lines for each of the 5 sources observed with NIFS. Left to right: The Ha intensity, Ha velocity and Ha 
velocity dispersion maps, and the available imaging (we show HST imaging if it is available otherwise we display the ground based K band imaging). For 
SMMJ0332— 2755 imaging in the i,j and k bands are available therefore a 3 colour image is displayed. The boxes marked by the red dashed lines represent 
the field of view of NIFS on the imaging. Top to bottom: RGJ0331-2739, SMMJ0217-0505, SMMJ0333-2745, SMMJ22 18+0021, SMMJ0332-2755. 
The green dotted lines mark the division into components, labelled as 'a' and 'b'. The yellow circles indicate the FWHM of the seeing dominated PSF from 
observations of the standard stars. This estimate from the standard stars is, however, likely better than what is achieved during the long science exposures. 
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Figure 4. Maps of the properties of the Ha emission lines for each of the 3 sources observed with SINFONI with strong Ho detections. Left to right: The Ha 
intensity, Ha velocity and Ha velocity dispersion maps, and the available imaging (we show HST imaging if it is available otherwise we display the ground 
based K band imaging or the IRAC 3.6fim channel map.). The boxes marked by the red dashed lines represent the field of view of SINFONI on the imaging. 
Top to bottom: SMMJ2217+0017, RGJ0332-2732 and SMMJ0217-0503. The green dotted line mark the division into components for SMMJ0217-0503. 
The yellow circles indicate the FWHM of the seeing dominated PSF from observations of the standard stars. This estimate from the standard stars is, however, 
likely better than what is achieved during the long science exposures. 



system rather than a merger, however this is dependent on the ori- 
entation, configuration and halo mass. 



4 RESULTS 

4.1 Dynamical Properties 

In order to test whether the SMGs are distinct in their basic dy- 
namical properties from other galaxy populations, we compare the 
kinematics of our SMGs to other high-redshift star-forming galax- 
ies studied in previous surveys. In particular, we calculate the ratio 
of the maximum velocity across the source to average velocity dis- 
persion, which provides an indicator of the extent to which a sys- 
tem is rotation or dispersion dominated, V b s /2cr. Here, V b s is 
the difference between the maximum and minimum velocities de- 
rived across the spatial extent of the galaxies, taken as the 5th and 
95th percentiles in the velocity distributions determined from the 



Ha line-fitting, while a is the integrated galaxy velocity dispersion 
measure derived from the integrated spectr a. In Fig. [5] we compare 
our res ults to the star-forming galaxies from lForster-Schreiber et all 
(2009), UV/optically selected galaxies in a similar redshift range to 
our sample of SMGs (the SINS sample). The SINS sample con- 
tains a wide variety of sources and hereafter comparing to this 
SINS sample refers to the core majority of the SINS sample of 
turbulent disk-like systems. We note that the integrated velocity 
dispersions will be larger than the average dispersions measured 
in the maps (Figs. [3] and @ since the integrated values contain a 
non-negligible contribution from the range of velocities within the 
sources. We compare the average velocity dispersion in the maps to 
the integrated value finding that, on average, the integrated values 
are ^20% larger than the average in the dispersion maps. However, 
since lForster-Schreiber et al.l ( 120091) use the integrated velocity dis- 
persion in their calculation of V b s /2cr we also use the integrated 
value to enable a consistent comparison. 
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Figure 5. Properties of the SMG sample compared to the sample of SINS galaxies. V bs/2cr is ratio of half the observed velocity gradient to the integrated 
velocity dispersion and provides a measure of how rotation or dispersion dominated a system is. We note that in the case of merging systems V b s /2cr is 
not a direct probe of the physical state of the galaxy since it does not represent the velocity gradient across a single galaxy, however it remains useful for the 
comparison of the kinematic fields to the SINS sample. The error bars represent the 1 <r uncertainties. Top left: Relation between the uncorrected Ha SFRs 
and V b s /2(T. Top right: Relation between the extinction corrected Ho SFRs for the SINS sample, the far-infrared SFRs for the SMGs and V b s /2<r. Bottom 
left: Relation between the Ha half light radii and V b s /2<T. Bottom right: Relation between the dynamical mass and V b s /2o\ The masses of the SMG 
sample represent a useful parameterization in terms of the integrated- velocity dispersion and rj m for comparison to the SINS sample, but only have a physical 
representation as a dynamical mass for certain model configurations in virialized systems. We divide SMMJ0217— 0503 into its two components and mark 
the northern SMMJ0217— 0503 source (SMMJ0217— 0503a) as an 'AGN' (open red diamond) and the southern disk-like component as a SMG. We plot the 
properties of the whole SMMJ0217-0503 system also marked as an 'AGN'. We highlight RGJ0332-2732 and SMMJ0332-2732 as 'candidate AGN' (open 
red squares) since there is some evidence that the velocity fields may be partially affected by AGN given their high values of [NII]/Ha (Section 3.3). 



Fig.[5]shows V bs/2<7 plotted against the SFRs, Ha half light 
radii and the dynamical masses. The average V bs/2<7 of the SMGs 
is 0.9±0.1 which is higher than the average of the SINS sample; 
(V obs /2o-) =0.60±0.05. However the high values of V obs /2a 
are due to large velocity gradients within the systems rather than 
especially low velocity dispersions. The average V b s of our 
SMGs is 330±80kms _1 compared to 160±20kms~ 1 for the SINS 
sample whereas the average integrated velocity dispersion of the 
SMGs (excluding the AGN, SMMJ0217-0503) is 130± 20kms _1 
which is consistent with the average for the SINS sample of 
130±7kms~ 1 . The large velocity gradients in the SMGs may be 
due to the minimum and maximum velocities coming from dif- 
ferent components of a merger system. Indeed, the systems which 
do not have clear sub-components in Fig|4](SMMJ2217+0017 and 
RGJ0332— 2732) have the lowest values of V b s /2<7, along with 
the broad-line AGN component SMMJ0217-0503a. The largest 
value of Vobs is in the southern component of SMMJ0217— 0503b 
suggesting a sharp velocity gradient through the disk or possibly 
a large velocity offset between merging components. When con- 
sidering merging systems, V bs/2cr does not represent the velocity 
gradient across a single galaxy, and is no longer a direct probe of 
the physical state of a turbulent disk. However, it remains a useful 



comparison of the kinematic fields of our SMGs to the SINS galax- 
ies as we attempt to ascertain how the SMGs differ from typical 
z ~ 2 disk galaxies. 

The average SFR H c« of the SMG sample (41±7M Q yr _1 ) is 
somewhat higher than (but within the range of) the SINS sample 
(17±2 Mq yr _1 ). Correcting the SMGs for the effects of extinction 
using A v = 2.9 ± 0.5, the ave rage A v calculated from a sample 
of SMGs bv lTakata et all [2006), raises the SFR values to approxi- 
mately align with the SFRfir values. In the top right-hand plot of 
Fig. [5] we compare the SFRfir values of the SMGs to the extinc- 
tion corrected SFRh q values of the SINS sample since we do not 
have individual extinction values for the SMGs. We find that the 
SMGs have higher SFRs, with (SFRfir) = 370 ± 90 M yr _1 , 
than the SINS sample which has (SFR Hq ) = 165 ± 35 M yr _1 . 
The value of A v used in the SINS sample correction comes from 
SED fitting and also takes into account the fact that Balmer emis- 
sion is more attenuated than starlight without a direct measurement 
of this extra attenuation therefore there are significant uncertainties 
in this correction factor. 

The Ha half-light radius measurements for the SMGs span 
a narrower range than the SINS galaxies, with (ri/2) = 3.7 ± 
0.8 kpc. The SINS sample have increasing V bs/2cr with ri/2 sug- 
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gesting that as the size of the galaxies increases the galaxies be- 
come more rotation dominated. Alternatively, as the sizes increase 
the velocity differences also increase perhaps due to the velocities 
having larger ranges when they span larger systems since they lie in 
more massive halos, however when considering merging systems 
the size is defined by the merger state rather than the halo mass. 
This offers an explanation for why we do not observe the strong 
relationship between Vobs/2cr and riy 2 in the SMG sample, since 
we interpret these systems to be mergers (Section [4.21 . 

The SMG Ha morphologies and velocity structures clearly 
show disturbed and complex dynamics, which motivates using a 
dynamical mass calculation assuming the systems are not dom- 
inated by rotation. We therefore use a virial mass approxima- 
tion, Equation [3] which is the method used to calculate dynamical 
masses f or the dispersion-dominated g alaxies of the SINS sample, 
stated bv lForster-Schreiber et al.l d 2009l) as being suitable for a wide 
range of galaxy mass distributions dBinnev & Tremaindl2008l) . 



M dyn = 



:7a 2 T 1 



/■I 



(3) 



Since we interpret the SMGs as being merger systems (Sec- 
tion |42]l and the velocity dispersions may be driven by the star 
formation or AGN (Section [5. 2\ , the virial assumption may not be 
valid. However, this 'mass' provides a useful comparison between 
the SMGs and SINS galaxies in terms of the integrated velocity 
dispersion-rx/2 parameter space covered. The masses determined 
from Equation|3]of the SMGs span a range from 1.4 x 10 10 M to 
1.6 x 10 11 Mq. This does not include the three AGN candidates, 
which have systematically larger line widths, since the masses cal- 
culated using this method may simply reflect the gas motions close 
to the AGN. Our non-AGN SMG sample shows a similar range to 
that found in the SINS galaxies, and thus the average velocity dis- 
persions of the SMGs do not obviously distinguish them as more 
massive, however the masses are further investigated in Section [5Jl 
The interpolated restframe //-band luminosities of our SMG sam- 
pl e suggest they have s imilar stellar masses to the SMGs studied 
in lHainline et alj l l201ll) (median ~7x 10 10 M Q ). This places the 
stellar masses of the SMGs in the upper third of the estimated stel- 
lar masses of the SINS sample which supports the possibility that 
these ULIRGs may be formed from SINS-type galaxies merging. 

To quantify the difference of the two samples in these quanti- 
ties we apply a Kolmogorov-Smirnov test (including only the non- 
AGN component of SMMJ0217-0503). We derive a 40% proba- 
bility that the SINS and SMG samples are drawn from the same 
distribution in terms of r 1 / 2 and 96% for Md yn suggesting that the 
SMGs are not significantly different to the SINS sample in terms 
of these properties. The probability that the SMGs' SFRpir are 
drawn from the same distribution as the dust-corrected SFRho in 
the SINS sample is 0.02%, while for V bs/2<7 the probability is 
2%, suggesting a difference in these quantities between our SMG 
sample and SINS galaxies. 

We also consider how the galaxies in the SINS sample classi- 
fied as mergers compare to our SMGs. The SINS merger subsam- 
ple are classified as mergers using kinemetry, the method we use in 
Section|4j2] Particularly in V bs/2a and SFR, the SMGS are dis- 
tinct from the SINS merger subsample, which may be indicative of 
the lower luminosity mergers found in the SINS galaxies. 

Despite the large difference in their SFRs, the SINS galaxies 
and SMGs can not be clearly differentiated simply from basic mea- 
sures of their velocity difference, velocity dispersion or half-light 
radius, and so more complex dynamical modelling is required. The 
SMGs could still represent a mixture of disks and mergers in the 



same way the SINS sample does, with SMGs being scaled up and 
more clum py versions of lo cal disks as seen in the simulations de- 
tailed bv lDaveetaf] ( l2010h . When we consider this further in Sec- 
tion l4.2l where the levels of asymmetry in the velocity and disper- 
sion fields are measured, we find strong evidence that the dynamics 
of all the SMGs are more consistent with being mergers. 



4.2 Evidence for Mergers and Kinemetry analysis 

The intensity maps of Fig. [3] and [4] show the distribution of star 
formation within the sources via the Ha morphologies. Six of the 
sources exhibit clear multiple peaks in Ha intensity, three of which 
match multiple peaks in imaging, which could in some cases be a 
result of structured dust within the sources obscurin g regions of star 
format ion resulting in clumpy Ha morphology 1 Swinbank et al.l 
l2010al) . The multiple peaks could also be distinct objects with are in 
the process of merging. The velocity maps show large ranges in the 
velocities of the systems (typically 400 kms -1 peak to peak) and 
the dispersion maps show disturbed motion within all the sources, 
providing initial evidence that most of the SMGs could be merger 
systems. Indeed only RGJ0332— 2732 (an AGN candidate) exhibits 
a smooth velocity and dispersion field. 

To better quantify the dynamics, we attempt to fit the 
two-dimensional velocity field using simple disk models. We 
parametrise the two-dimensional velocity field using an arctan 
function of the form V(r) = ^ V c arctan(r/r t ) where V c is the 
asymptotic circular velo city and rt is th e radius at which the ro- 
tation curve turns over jCourteaulll997l) . The model has six free 
parameters, disk centre [x/y], peak rotational velocity, V c , r t , disk 
inclination, ((') and sky position angle. To find the best fit, we use a 
genetic algorithm with 10,000 possible random fits in each gener- 
ation, and calculate the \ of eac h fit. Each parameter is then ad- 
justed towards the best fit parameter and the process is repeated un- 
til all the values have converged (and we demand >10 generations 
have passed before testing for convergence). Unsurprisingly, given 
the complex dynamics seen in Figs.[3]and[4] simple disk models do 
not provide an adequate description of the data, with average dif- 
ferences between the best-fit model and the data of 40±10kms _1 
and a median reduced % 2 of 60±30. As a control sample, we also 
run the same fitting procedure on a sample of disks at high-redshift 
(Swinbank et al. in prep.) finding a median rms of 27±7 kms~ 1 and 
reduced \ 2 values between 5 and 8 which suggests the SMGs are 
not consistent with being disk-like systems. 

In order to provide a more definitive test of whether the kine- 
matics are rotation versus merger dominated, we can use kinemetry, 
which was originally developed to search for non-circ ular dynami- 
cal fe atures in the stellar dynamics of local ellipticals (Co pin et al.l 
l200ll) . but has been extended to t he dynamics of high- z galaxies. 

The kinemetry described in Sha piro et al.l ((2008) (employing 
the code ht tp://www.e so.0rg/~dkrajnov/idl/l uses the method de- 
tailed in I Krai no vie et alj ( I2006J). and has been used to ident ify the 
mergers in the SINS sample jForster-Schreiber et al.l 12009). This 
might allow us to better test whether the dynamics of our SMGs 
are different from SINS galaxies. The asymmetry is measured by 
fitting ellipses, varying the position angle and ellipticity, to the ve- 
locity and dispersion fields. The ellipses are fit at increasing radii 
and the coefficients of the Fourier expansion of each best fit ellipse 
at each radius are recorded. These coefficients are then used to es- 
tablish the level of asymmetry level. The velocity field asymmetry 
is measured as the average over all radii of 
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Vasym — ( t-> ) r (4) 

£>i,v 

where k avg , v = (k 2 , v + k 3iV + k 4 , v + ks, v )/4 and k n , v = 
v/ An, v + B„,v and where A lljV and B njV are the coefficients of 
the Fourier expansion of the best fit ellipse. The velocity dispersion 
field asymmetry is measured as the average over all radii of 

_ / kavg.j , ... 
C'asym — \ -r> / r y J f 
Dl,v 

where k avgiCT = (ki >CT + k 2 , CT + k s , CT + k 4 , CT + k 5 , CT )/5 and 
k Ilj(T = y An, a- + Bi jCT and where A n . CT and B nj(T are the coeffi- 
cients of the Fourier expansion of the best fit ellipse. 

The two asymmetry measures are combined together as 
K asym = \/v asym + aisym. Systems with K asym > 0.5 are 
classified as mergers. This limit is determined by running the 
ki nemetry analys i s on template merger and disk systems, detailed 
in Ishapiro et al.l | |2008|) . to find the two populations separate at 

K a sym ^ 0.5. 

We fit the centres of the ellipses first and then we re-run the 
kinemetry 100 times per source varying the centres within the stan- 
dard deviation of the initial centre fit. The v asym and <7 asym val- 
ues we quote in Fig. [6] are the bootstrapped median and standard 
deviation of these 100 runs. The K asym values we quote (in Ta- 
ble^ are the bootstrapped median and standard deviation of the 
100 runs (rather than using the median v asym and o- asym values). 
Owing to the lower quality of the velocity and dispersion maps of 
SMMJ0332— 2755 we bin the maps by 2x 2 pixels in order the run 
the kinemetry. 

In Fig. [6] we compare the levels of asymmetry in the veloc- 
ity and dispersion fields to a number o f observed and simulated 
systems, detailed in Shapi roet al.l d2008h . The SMGs have, on av- 
erage, higher levels asymmetry that both the low-redshift ULIRG 
population and the more 'normal' star forming galaxies of the SINS 
sample. The SMGs also do not lie close to the SINGS low-redshift 
spiral galaxies confirming that the SMGs are very different to other 
populations of star forming galaxies in terms of their gas dynamics. 

The background red-blue pixels represent the results of run- 
ning the kinemetry o n the template disks and mergers provided by 
Shapir oet al.l d2008l) . We find all of the systems in our sample of 
SMGs have K asym > 0.5 (Tabled and thus we classify all systems 
as mergers. We also calculate the level of asymmetry in the south- 
ern component of SMMJ0217-0503 finding that K asym < 0.5 
which suggests that this component is a disk-like source (as im- 
plied in Figs. [2] and |4), merging with the northern AGN source. 
Indeed the K asym value for the entire system is 3.9±0.3 and there- 
fore the whole system is clearly classed as a merger. The sources 
which are the closest to the disk region of Fig.|6]are the AGN can- 
didate source (SMMJ2217+0017), RGJ0331-2739 and the disk 
component of SMMJ0217— 0503 (the southern component in a ma- 
jor merger). However the Kasym values for SMMJ22 17+00 17 and 
RGJ0331— 2739 (using the bootstrapped median of the 100 runs) 
put these sources in the merger class. Therefore we classify all sys- 
tems as mergers. 

The velocity and dispersion maps cover, on average, a total 
area of ~160kpc 2 (at z — 2). The spatial resolution of the obser- 
vations is 0.5" corresponding to an area of ~16kpc 2 therefore on 
average we have ~ 10 resolution elements in each map. We perform 
an identical analysis on our SMGs to that done on the SINS sam- 
ple, which have similar spatial extent. We note however that the 
fit is formally under-constrained and as such, where we find a disk 
solution, we cannot completely rule out an underlying merger. 



We also plot the asymmetry values determined from running 
the kinemetry on a sample of simulated SMGs at z = 2 from 
iDave et all (2010). We find these simulated SMGs have a range of 
asymmetry values forming a sample of both disks and mergers and 
our observations are not in complete agreement with the simulated 
systems since we do not observe the range of disks predicted by the 
simulations. We perform a two dimensional KS test on the model 
SMGs and observed SMGs to find the probability that both samples 
are drawn from the same distribution of 0.1 and therefore we can- 
not rule out that the simulated SMGs represent the observed SMGs. 
Further, the merger-like kinematics in these simulations does not 
necessarily imply that the high infrared luminosity is being driven 
by a merger event, since some of these simulated galaxies are not 
actively merging. 

There is therefore considerable evidence for the SMGs to be 
merging systems; the multiple peaks in the Ha intensity fields, the 
disturbed velocity and dispersion fields to which 2d disk models 
could not be fit, the distinction in some integrated quantities be- 
tween the SMGs and the SINS star forming galaxies, and finally 
the kinemetry analysis which classes the SMGs as mergers. 

4.3 Sub-components and separations in the Ha intensity field 

Having established that the Ha dynamics of the SMGs are likely 
merger-driven, we next try to identify the components which are in 
the process of merging, as we have done already for the obviously 
distinct SMMJ0217-0503 a & b. We noted before that some of 
the systems have clearly separated Ha intensity peaks, and in more 
than half the sample these clearly align with structured velocity dis- 
persion peaks, as would be expected for merging, self-gravitating 
components. The sources are therefore split into components by 
identifying the minima in the intensity maps, placing the two in- 
tensity peaks in different components as shown by the green dotted 
lines in Figs. [3] and [4] 

The components of SMMJ22 18+0021 and SMMJ0332-2755 
also match the components seen in the HST imaging shown in Fig. 
[3] We find the division is possible in all sources except the AGN 
candidates, SMMJ22 17+00 17 and RGJ0332-2732. It should be 
noted that splitting the sources into components does not necessar- 
ily represent two distinct objects which are merging since these sys- 
tems could be in the later stages of merging where multiple passes 
have mixed the gas and stars. However, we discuss the component 
properties as a reference point assuming that we have identified 
real merging components, acknowledging that these systems may 
reflect more complex dynamics than a simple merging pair. 

A two-dimensional disk model is well fit to the southern 
component of SMMJ0217-0503(b) with a reduced \ 2 of 1.4 and 
higher order moments (kinemetry) also clearly identify this south- 
ern component as a coherent disk (detailed in Section 14.2) . How- 
ever we also note that two regions within SMMJ0217— 0503b (the 
eastern and western halves) show different (low) ratios of [NII]/Ha 
(0.44±0.07 and 0.25±0.07 for the eastern and western sections re- 
spectively). If not dominated by ionization effects, then this sug- 
gests that this component may contain regions of varying metallic- 
ity and therefore could itself be the product of a recent merger. 

We derive V bs/2<r, I1/2, Md yn and SFRh q for the compo- 
nents in the same way as for the whole systems (shown on Fig.|5j- 
The components exhibit similar V bs/2cr values to the entire merg- 
ing systems since Vobs and a are both smaller in the components 
((V oba )=260±40kms" 1 and (cr)=110± 10 km s _1 ). The compar- 
ison of our SMGs' merging components with SINS galaxies sug- 
gests it is possible that the merger of two typical SINS isolated disk 
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Figure 6. The asymmetry measures of the velocity and velocity dispersion fields for the sample of SMGs compared to a number of other samples; the SINGS 
sampl e of local spiral galaxies, a sample of low-redshift ULIRGs, the SINS sample of z ~2 star forming galaxies, detail s of which are given bv lShapiro et alj 
1 2008 ). The background red-blue points are the results of running the analysis on template disks and mergers taken from Shapiro et al. I 2008j). All of the SMGs 
lie in the merger region (red background) classifying all of SMGs as mergers. We highlight the AGN candidates, RGJ0332— 2732 and SMMJ22 17+00 17 and 
the whole SMMJ0217— 0503 system containing the AGN and disk and we also plot the asymmetry values of SMMJ0217— 0503b only (disk-like component 
of the merging system). We plot the results of running the kinemetry on a sample of simulated SMGs at z = 2 finding they represent a range of disks and 
mergers. 



galaxies may generate a SMG. However, the properties of the in- 
dividual components will develop throughout the merger, evolving 
with each pass and therefore our division into components may not 
correctly separate the two progenitor merging sources. 

In Fig.|7]we show the spatial offsets for our sample, as well as 
for a sample of SMGs obser ved in CO, and prev ious Ha studies. 
These are compared to local dVeilleux et al.l l2002) ULIRGs which 
have multiple components with small separations (<4kpc apart), a 
few having components up to 25kpc a way. This matches merger 
simulations (Mi hos & Hernquisti 1996) where the large burst of 
star formation occurs at the final stage of the merger and thus the 



ULIRGS are expected to be observed in very close configurations. 
This is however not observed in the SMGs with most of the Ha 
sample having separations greater than 4kpc indicating that the 
ULIRG phase can occur at the early stages of merging in these 
high-redshift systems. The large separatio ns observed in CO are 
predicted by hydrodynamical simulations IIHavward et al 1 l2011al) 
which explain the 850 /im source population. 

These 3 samples are complete with the local ULIRGs being 
cut at Lir = 10 12 L Q and both the CO and Ha SMG samples 
spanning a narrow range of Lir = 2 — 10 x 10 12 Lq. There is 
however a selection bias since sub-components with r<4 kpc would 
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Figure 7. Histograms showing the distributions of the projected separa- 
tions of the SMGs (Ha observations in the top plot and CO observa- 
tions in the middle plot) and a sample of local ULIRGs (bottom). The 
Ha observations combine the NIFS and SINFONI samples with all ex- 
isting IFU s ample s detailed in Swinbank et al. 12006) and also including 
iTecza et al.l j2004l) . [Swinbank et alj J2005l) and lMenendez-Delmestre et alj 
<2012t) . RGJ0332— 2758 from Casey et al. (in prep.) is also included. The 
Ha studied SMGs have larger separations on average than local ULIRGs 
suggesting different merger stages highlighted by each population. This 
large separation is not seen in the CO observed SMG sample potentially 
highlighting the differences in the gas distributions traced by Ha and CO. 



not be resolved in the Ha data. Furthermore, the NIFS Ha sample 
is not sensitive to separations greater than 15 kpc as the field of 
view of NIFS is approximately 30 kpc by 30 kpc (at z — 2) whereas 
the CO SMG sample could detect separations up to ~30kpc. The 
maximum separation detectable by SINFONI is 30 kpc at z — 2 
due to the larger total field of view of SINFONI, however we nod 
the source around 4 quadrants of the IFU and therefore the effective 
field of view is approximately 30 kpc by 30 kpc. 



5 DISCUSSION 

The multiple Ha peaks, the failed disk-fits to the velocity fields and 
the kinemetry analysis together provide considerable evidence that 
these systems are likely to be major mergers. We directly identify 
multiple components in the majority of the systems and suggest that 
these are in the process of merging. Our findings agree with tho se of 
ISwinbank et alj J200r3) and lMenendez-Delmestre et who 
also find multiple-component complex systems in their Ha obser- 
va tions of SMGs. Th ese findings align with the conclusions made 
bv lEngel et alj < f20ich . whose CO obse rvations suggest that SMGs 
are often formed in m ergers, and also iRiechers et al.l (1201 lab and 
Riec hers et ai]( l2011bl) which detail EVLA observations of 3 SMG 
systems in CO(J=l— 0) presented as either early stage or late stage 
mergers. Next we further explore the multi-component nature of 
the SMG systems and study the observed offsets in the context of 
the parent dark matter halos. 



5.1 Masses from component properties 

SMGs have a clus tering length r ~ 8Mpc dBlain et all 1 20041 ; 
Iffickox et alj l2012h consistent with a form of evolution ensur- 
ing their properties subsequently match the clustering of z ~ 1 



evolved red galaxies, and finally of galaxy clusters near the present. 
This should indicate that they occupy the largest dark-matter halos 
(~ 10 13 M Q ). In order to examine whether the multiple compo- 
nents of the SMGs can be used as a test of the halo mass we inves- 
tigate the relationship between the projected separations and the 
velocity offset of the components. 

Fig.[8]shows the spatial and velocity offsets of the components 
within each of our systems, combined with thos e from previous 
SMG Ha-IFU and CO observations summarized in Swinbank et al. 
(2006), a long with two additional SMGs w ith Ha-IFU measure- 
ments in iMenendez-Delmestre et alj d2012l) and Casey et al. (in 
prep.). We also plot the two SMGs in our sample which cannot be 
divided into components (the AGN candidates SMMJ22 17+00 17 
and RGJ0332-2732) with dV, dR limits. The velocity offset for 
these sources is the difference between the minimum and maxi- 
mum velocities in the systems (Vobs) and the spatial offset is ri/ 2 , 
since we assume that if there was a recent merger, it has progressed 
so that the initial systems are less than rx/2 apart and the velocity 
difference must be below Vobs- 

Fig. [8] highlights that the CO measurements have generally 
found significantly smaller component separations than this com- 
bined sample of IFU Ha maps ((dR) = 4 ± 1 kpc for the CO sam- 
ple compared to (dR) = 10 ± 3 kpc for the complete Ha sample), 
while the CO velocity offsets extend to high values, with an average 
of 440±90 kms -1 compared to 230±40 km s" 1 for the Ha s ample , 
consistent with virialised halo models from Swin bank et al. (2006). 
These clear differences between the CO and Ha samples may be 
indicative of the sensitivity of Ha to lower levels of star formation 
offs et from the bright CO region. Ind eed, ISwinbank etail d2005l) 
and IMenendez-Delmestre et alj ({2012) have clearly demonstrated 
this to be the case for these two SMGs with both measurements. 

Both the spatial and velocity offsets found from the Ha com- 
ponents in our SMGs are similar to that of the other SMGs observed 
in Ha with IFUs ((dR) = 7±lkpc, (dV) = 180±50kms _1 and 
(dR) = 8±2kpc, (dV) = 200±100kms- 1 for the literature sam- 
ple and the sample in this work respectively). SMMJ0217— 0503 is 
an outlier in the plot with its extremely large spatial and velocity 
offsets. Given the similarities in overall properties of these vari- 
ous SMGs, we cannot attribute this difference to anything other 
than small sample variance. However, the wider range of dV found 
in our expanded SM G sample does make a simple virialised halo 
model, as adopted in lSwinbank et al. a less satisfactory de- 

scription of this dataset. 

Our goal is to estimate dynamical halo masses for these galax- 
ies, but this is non trivial for what then may often be unvirialised 
systems, and so we turn to N-body simulations and extract the av- 
erage dR a nd dV for a set of ma ssive halos from the Millennium 
simulation (Sprin gel et alll20 05b). In Fig.[8]we show the projected 
spatial and velocity offsets for bound sub-halos within dark mat- 
ter halos of total mass (M). All halos within each mass range are 
identified and the position and velocity of the most bound sub-halo 
within each is extracted. All remaining sub-halos are then identified 
and the spatial (dR) and velocity (dV) offsets between these sub- 
halos and the most bound sub-halo are calculated for each mass 
range. We then bin the data in bins of dV and plot the median in 
each bin to create the tracks for each halo mass. 

The mass range of 13 < log(M[ii -1 M0]) < 14 appears 
to best represent the combined dataset. For comparison, we show 
the interpolated relationship to log(Mhaiofh~ *M fl l) = 12 8t°'l , 
the SMG halo mass measured from clustering dHickox et al.l20120 . 
which is lower than the average halo mass we derive. The clus- 
tering measurements however come from relatively small sam- 
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Figure 8. Relationship between the offset in position and velocity of the components within the SMGs. Included are the limits on the separations for the 
two sources which cannot be divided into components (the AGN candidates SMMJ22 17+00 17 and RGJ 0332— 2732). The outli er is the AGN-disk system of 
SMMJ0217— 0503 with a large velocity offset between the two systems. We include the offsets detailed in Swinbank et al. (2006) fro m the i r IFU observations, 
and also their compilation of other IFU offsets and those measured from lon gslit and CO observations llFraver et al.ll998tlNeri et alj2 003; Genz eFet alj2003l ; 
iKneib et alJl2004l : lGreve et al^OOjljTacconi et all200fj) . SMMJ123549.44 llvIenendez-Delmestre et alj2012h and RGJ0332-2758 (Casey et al. in prep.) are 
also plotted. The dashed lines represent the median projected offset between the most-bound sub-halo and the other subhalos within dark matter halos in the 
quoted mass range from the Millennium Simulation (Springel et al. 2005b). The dotted line represents the interpolation to log(Mhalo[h — 1 Mq]) = 12.8^Q 5, 
the SMG halo mass measured from clustering (Hickox et al. 2012). The observed offsets and the measured offsets from the simulation database are not 
corrected for viewing orientation of the randomly orientated systems. 



pies of SMGs and therefore have large statistical errors. Further- 
more, the bias of SM Gs to halos is potentially quite complicated 
( Chap man et aTT2 0Q9). given the short duration of the burst and the 
unknown range of environments that they occur in. 



The large spatial and velocity offsets of SMMJ0217-0503 
suggest this system could have a log(M[h _1 M0]) ~ 14 halo 
which we view along the line of sight velocity, or an even larger 
halo mass if it is not in this preferred viewing orientation, providing 
evidence that SMGs may sometimes be beacons fo r massive proto- 
clusters of galaxies (loR(M[h~ 1 Mp ) 1) > 14) dChapman et al.l 
120091 : iTamura et all l2009t iMatsuda et all l201lh . In the Millen- 
nium Simulation database there are 28 halos with masses 14 < 
log(M[h. M0]) < 15 confirming that these large masses are pos- 
sible and therefore that the two components of SMMJ0217— 0503 
are likely part of a common halo rather than two separate sources 
whose projection places them within our field of view. However we 
note that it is also possible that the SMMJ0217— 0503 system may 
represent a flyby rather than a merger owing to the large velocity 
offset of 670±70kms _1 . 



5.2 High velocity dispersions and the relationship to SFR 
density 

Figs.[3]and|4]clearly show resolved substructure in the Ha velocity 
dispersion maps which often correlate with the intensity fields. We 
further investigate this relationship and explore the origin of the 
high velocity dispersions found in our SMGs. Despite the larger 
V bs/2cr in the SMGs compared to the SINS galaxies, the values 
are significantly lo wer than observed in local disks (v/cr ~ 2 — 5: 
Hunter et alj2005h which have much lower dispersions but compa- 
rably large circular velocities. Given our findings (in Section [42} 
that the asymmetrical velocity and dispersion fields enable us to 
classify the systems as mergers, it is natural that the high velocity 
dispersions within our sample may be due to the merging systems, 
causing the gas to be disturbed. 

The SINS galaxie s with high velocity dispersions are explored 
in iGenzel et alj d2008h within the context of disk-like galaxies dis- 
playing turbulent motion. They explain the turbulent, rotating disk 
structures either usin g feedback from the star formation driven by 
supernovae winds (Efstat hioulboOOl) . or with accreting gas driv- 
ing the turbulence as it flows into the disk from the dark halo 
(Genzeletal. 200^). By contrast, the simulations of iDave et al.l 
i 201Ch suggest turbulent motion in massive disk-like galaxies is 
driven by harassment by smaller satellites and the infall of gas. In 
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the SMGs, the high Ha velocity dispersions may also be connected 
to the star formation itself, rather than the forces acting on the gas 
from the merger. 

We first test whether the relationship between the velocity dis- 
persions and the star formation rates can be explained as being 
purely a result of the Kennicutt-S chmidt (KS) star formation den- 
sity - gas relation (Kennicutt 1998b): 

Ssfr = ASg as (6) 

We relate E gas to the velocity dispersion (a) by setting 
Toomre's Q Criterion to equal to 1, since we assume the gas finds a 
stable disk configurati on in a short time as suggested by hydrody- 
namical simulations ( Springel et al. 2005a), where 



and k = \/2V/r dSchave& Aguirrd 120051) where we take 
V = V obB and r = ri /2 . We ta ke A to be 10" 3 - 63 and n = 1.27 as 
derived in lGenzeletal.ld2010l) . 

We plot this relation in Fig. [9] to enable a comparison to the 
observed relations between the Ha velocity dispersion per pixel 
and the star formation intensity. We note first that there are several 
selection effects at work. There are sections of the a — Esfr plane 
which cannot be populated from the observations. The upper left 
region of Fig.|9]is empty as broad lines with low star formation in- 
tensities have too low S/N for Gaussian fitting. There is also a lower 
limit to the velocity dispersion restricted by the spectral resolution 
and the lower region (~< SOkms" 1 ) cannot be populated. 

The high velocity dispersion per pixel values of our three 
SMGs showing some AGN characteristics could be arising from 
ionization of gas in the broad- and narrow-line regions of the AGN 
itself and it is possible that the lower Ha densities in these sources 
are due to feedback from the AGN i nhibiting the star form ation 
in the regions surrounding the AGN (Nesvadba et al. 2008). One 
SF-dominated SMG, SMMJ0217-0503b, which has intrinsically 
narrow lines, is not sufficiently resolved spatially to separate the 
double-peaked rotation profile in the central region, and our line- 
fitting procedure blends this double peak as a single line. We there- 
fore choose three distinct regions within the source and use the in- 
tegrated spectra of these regions to obtain the velocity dispersion 
and star formation rates, allowing a fit to the double peak. 

Fig. [9] shows that the KS relation does not describe the data 
well since the majority of the data have higher values of velocity 
dispersions at a given star formation intensity than predicted by 
this model. Clearly, both the dust extinction and inclination will 
have a dramatic effect on the location of these tracks. An average 
A„=2.4±0.6 is required for the SMGs (excluding the AGN candi- 
dates) to match this Esfr-c relation, which is compati ble with the 
average implied dust correction derived for SM Gs (eg ISmail et al.l 
l2004| - |Swinbank et alj2004l ; lTakata et al.l2006h . and suggestive that 
SMGs simply follow the s tandard KS-law (as suggested may be the 
case bv llvison et alj201lh . 

Since we cannot constrain the dust extinction or inclination 
angle from our data we also consider if the velocity dispersions are 
driven by the star formation i tself. We follow the same analysis as 
carried out in iLehnert et al.l (120091) when exploring the source of 
the turbulence in SINS galaxies. The SMGs are comparable to the 
SINS galaxies, with the majority of the data l y ing w ithin the range 
spanned by the relations from ILehnert et al. (2009) which repre- 
sent the plausible ranges in efficiency of the energy transfer from 
the star formation to the gas. It is therefore possible that it is the 
star formation which drives the velocity dispersion. However it is 



also possible that the turbulence is driven by the merging processes 
which stir up the gas. We require extinction maps to establish the 
true star formation intensities in order to further test these hypothe- 
ses. 



6 CONCLUSIONS 

We analyse subarcsecond resolution Ha line data, taken using the 
integral field spectrometers Gemini-NIFS and VLT-SINFONI, of 
nine SMGs with 2.0 <z< 2.7. The gas dynamics and morpholo- 
gies were mapped by tracing the Ha line within the eight of the 
sources. We compare the velocity fields in this sample to a sample 
of UV/optically selected star forming galaxies at the same redshift 
(studied in the SINS survey) finding that the SMGs are kinemati- 
cally distinct, with higher SFRs and larger velocity gradients across 
the systems. The majority of the SMG sample are not consistent 
with being disk-like systems with multiple peaks in Ha intensity 
and irregular and turbulent velocity and velocity dispersion fields. 
This is confirmed from the kinemetry analysis which classes all 
systems as mergers rather than disks. We identify multiple compo- 
nents in six of the sources further confirming the sug gestion that 
most SMG star bursts are triggered in major mergers (En gel et"ai1 

l2oioh . 

We bring together all existing IFU Ha observations of SMGs 
and use the spatial and velocity offsets between the components 
of the merging systems to infer an average halo mass, 13 < 
log(M[h _1 MQ]) < 14, for the systems which is higher than the 
mass predicted from clustering measurements. 

We observe higher values of V b s /2cr in the SMG sample 
compared to the SINS galaxies due to the large V bs values mea- 
sured across merging components. The velocity dispersion values 
in the SMGs are larger than observed in local disks. We explore 
the origin of this turbulence finding that the disturbed gas motions 
are consistent with the Kennicutt-Schmidt law and variable extinc- 
tions, but might also be driven by the merging torques or even the 
star formation itself. 

Analysing the outputs from the various simulation of SMGs 
in an identical way to the observations will enable us to better de- 
termine which situation best describes the trigger for the ULIRG 
phases in these galaxies. Initial evidence from the kinemetry analy- 
sis of the simulated SMGs suggests that the simulated systems are 
a mixture of disks and mergers. By better understanding the trig- 
gers of the ULIRG phases we can also probe the build of mass in 
some of the most massive galaxies in the Universe and investigate 
if these SMGs ar e the progenitors of the massive elliptical galaxies 
observed today (Hainline e t alj|201 Tl i lGonzalez et alj201lh . 
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Figure 9. Velocity dispersion against the star formation intensity per pixel for each source. The dot-dashed line represents the relation for the case i f the velocity 
disper sion values are purely a result of the Kennicutt-Schmidt (KS) relation. The dashed, dotted and solid lines represent models, detailed in iLehnert et al J 
1 2009), describing the relation between the velocity dispersion and the star formation intensity per pixel if the energy from the star formation is powering the 
motion of the gas. The three lines assume three efficiencies of the coupling between the input energy by stars and the ISM. The sources follow tracks along 
the plot despite being offset from each other. The SFRs are not extinction corrected. All sources are smoothed the match the seeing of 0.6" and then binned 
so that each pixel is 0.25" by 0.25" except for SMMJ0217-0503b. Since we observe the disk-like double-peaked Ha profile in SMMJ0217-0503b we plot 
the velocity dispersion and star formation intensity in three distinct regions across the source instead of using the results from the line-fitting procedure which 
blend the double peak as a single line 



HerMES is a Herschel Key Programme utilising Guaranteed Time 
from the SPIRE instrument team, ESAC s cientists and a mis - 
sion scientist. HerMES will be described in lOliver et"aL I d2012t) . 
The HerMES data was accessed through the HeDaM database 
(http://hedam.oamp.fr) operated by CeSAM and hosted by the Lab- 
oratoire d'Astrophysique de Marseille. 
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